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Abstract: Calsequestrin (CASQ) is a major Ca^^-storage/buffer protein present in the 
sarcoplasmic reticulum of both skeletal (CASQl) and cardiac (CASQ2) muscles. CASQ 
has significant affinity for a number of pharmaceutical drugs with known muscular 
toxicities. Our approach, with in silico molecular docking, single crystal X-ray diffraction, 
and isothermal titration calorimetry (ITC), identified three distinct binding pockets on the 
surface of CASQ2, which overlap with 2-methyl-2,4-pentanediol (MPD) binding sites 
observed in the crystal structure. Those three receptor sites based on canine CASQl crystal 
structure gave a high correlation {R = 0.80) to our ITC data. Daunomycin, doxorubicin, 
thioridazine, and trifluoperazine showed strong affinity to the SI site, which is a central 
cavity formed between three domains of CASQ2. Some of the moderate-affinity drugs and 
some high-affinity drugs like amlodipine and verapamil displayed their binding into S2 
sites, which are the thioredoxin-like fold present in each CASQ domain. Docking 
predictions combined with dissociation constants imply that presence of large aromatic 
cores and less flexible functional groups determines the strength of binding affinity to 
CASQ. In addition, the predicted binding pockets for both caffeine and epigallocatechin 
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overlapped with the SI and S2 sites, suggesting competitive inhibition by these natural 
compounds as a plausible explanation for their antagonistic effects on cardiotoxic side effects. 

Keywords: CASQ; calsequestrin; SR; sarcoplasmic reticulum; CASQ2; cardiac 
calsequestrin; CASQl; skeletal calsequestrin; TFP; trifluoperazine; ITC; isothermal 
titration calorimetry; daunomycin; DAN; desipramine; DSP; diltiazem; DTZ; propafenone; 
PFN; thioridazine; TDZ 



1. Introduction 

Clinical usefulness of many effective drugs is often limited by their toxicity as a consequence of 
dose- and time-dependent accumulation within muscle tissue [1]. A wide spectrum of possible 
mechanisms of those drug-induced muscle toxicities have been proposed, including muscle enzyme 
defects; changes in mitochondrial function and intracellular metabolism; the influence on membrane 
stability; and interactions with P-glycoprotein or Cyt P450 [2,3]. However, the precise 
pathophysiologic and biochemical mechanisms of most drug-induced muscle toxicities remain unclear. 
The main SR proteins responsible for Ca release process of the excitation-contraction coupling in 
cardiac and skeletal muscle are the Ca^^ release channel (RyR), junctin, triadin and the 
Ca^^-storage/buffer protein, calsequestrin (CASQ) [4]. There are several reports on the affinities of 
various drugs to RyR and their effects upon binding, such as dantrolene and caffeine [5-7]. In addition, 
both the skeletal and cardiac isoforms of CASQ (CASQl and CASQ2, respectively) show significant 
affinity for several classes of compounds such as anthracyclines, phenothiazines, and tricyclic 
antidepressants [8-10]. Many of those drugs can deposit in the sarcoplasmic reticulum (SR) and 
significantly perturb its critical role of storage and/or release of Ca^^ ions [10,11]. The structures of 
both CASQl and CASQ2 are composed of three thioredoxin-like domains [12-15], which, 
coincidently, is a binding motif for small molecules in many proteins [16,17]. Considering the high 
concentration of CASQ in muscle tissue (-100-500 mg mL"^), long-term or high-dose administration 
of those drugs could lead to their accumulation in the SR, despite their moderate affinity for CASQ 
(1-50 |LiM), which could substantially affect the Ca content of the SR due to decreased 
buffering/storage capacity of CASQ brought upon by its binding of these drugs [9,10]. In particular, 
those compounds or their metabolites could affect patients who already suffer from compromised drug 
metabolism and clearance or from impairments in proper CASQ function, such as those afflicted with 
catecholaminergic polymorphic ventricular tachycardia (CPVT) or malignant hyperthermia (MH) [9,18]. 

As has been noticed for flavonoids [19] and PDE-5 inhibitors [20], caffeine has been known to 
improve the survival rate of cancer patients, probably by modulating the toxic side effects of 
anthracycline-derivatives like daunomycin (DAN) and doxorubicin (DOX) [21]. In addition, 
epigallocatechin (EGC) and epigallocatechin-3-gallate (EGCG) have been reported to reverse 
DOX-induced intra-luminal Ca^^-depletion in the SR [22]. In this study, we applied molecular 
docking, isothermal calorimetry, and crystallographic techniques to find the potential binding sites. 
The structural mechanism behind the apparent affinities of CASQ for several drug molecules and a 
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structural model that explains the potential antagonistic effect of CAF and EGC were investigated 
through those comprehensive approaches. 

2. Results and Discussion 

2.7. Development of a Highly-Predictive Molecular Docking Model for CASQl-Ligand Binding 

We created an initial receptor model using the crystal structure of human CASQ2 (hCASQ2, PDB 
ID: 2VAF) to correlate the docking scores of drug compounds with their reported values from 
previous ITC studies [9]. For this purpose, we used a blind docking approach by including the whole 
protein in the affinity grid because the location of ligand-binding pockets in CASQ is unknown. The 
studied drugs were grouped into the following three binding classes based on their ITC dissociation 
constants (^d) (Figure 1): strong, moderate, and poor. The strong binders were amlodipine (AMD), 
daunomycin (DAN), doxorubicin (DOX), thioridazine (TDZ), trifluoperazine (TFP), and verapamil 
(VPM) the moderate binders were amitriptyline (AMT), diltiazem (DTZ), imipramine (IMP), 
nortriptyline (NTP), and promethazine (PMT); and the poor binders were chlorpromazine (CPM), 
ephedrine (EPH), isoproterenol (IPT), metoprolol (MTP), miconazole (MNZ), and propafenone (PFN). 
The competitive inhibitors, caffeine (CAF), epigallocatechin (EGC), and epigallocatechin gallate 
(EGCG), were also studied along with the known hCASQ2 binding drugs. 

Figure 1. The 2-D line drawings of the strong binders for human cardiac muscles 
(hCASQ2) and their corresponding experimental dissociation constants (^d)- Chemical 
structures for the complete set of the tested strong, intermediate, and poor binders are 
included in the Supplemental file 1 . 




Daunorubicin Doxorubicin Trifluoperazine 

Kd = 20.3|aM Kd=14.4|LiM Kd=18.2nM 




Amlodipine Verapamil 

Kd = 20.4|iM Kd= 14.8|LiM 



But this in silico model yielded a poor correlation between docking scores and experimental ^d's 
{R = 0.07), most likely due to potential inaccuracies in the side-chain conformations (Figure 2A) 
inherent to the low-resolution (3.8 A) of the protein structure. The same docking procedure was then 
repeated using the 2.4 A (PDB ID: ISJI) resolution crystal structure of canine CASQ2 (cCASQ2) in 
order to mitigate the effects of uncertainties in those side-chain conformations/orientations. This 
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significantly improved the correlation for drugs with strong or moderate affinities for hCASQ2 
(/?^ = 0.8) (Figure 2B). However, correlating the scores of poor binders with their corresponding 
experimental ^d's remained a challenge, as many of their scores were elevated to the moderate binding 
range, with some elevated to that of strong binders. Although the poor binders could not be correlated 
to the experiment, they were predicted to bind in the same pockets as those of strong and moderate 
binders. Visual analysis revealed less possibility for the presence of unique binding pockets for poor 
binders. This suggests that the bias in predictions might be coming from the scoring function instead, 
which could be over-predicting the binding interactions of poor binders. Inspection of the docking 
results for these poor binders highlighted that their aromatic cores were significantly smaller than those 
of moderate and strong binders, and possessed long highly flexible side chains with charges 
(Supplemental file 1). Hence, the large aromatic cores, flexibility of the side-chains and their sizes 
could be the crucial chemical features required for a tight fit into the CASQ2 binding pockets. 
Therefore, we excluded the docking predictions for poor binders on assumption that they suffer some 
docking bias and only focused on the strong and moderately binding molecules that could be well 
correlated to experimental values. 



Figure 2. The 2-D correlation plots of the molecular docking scores (AGb) versus the 
logarithm of experimental dissociation constants. Experimental log(^d) through 
least-square fitting using the crystal structures of hCASQ2 (2VAF-left, A) and cCASQ2 
(ISJI-right, B) as receptor models. 
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2.2. Prediction of Small Molecule Binding Pockets in CASQ2 



2_|_ 

Although binding of several pharmaceutical drugs and their side effects on both Ca -binding 
capacity and polymerization of CASQ have been well established [9,10], their binding sites and the 
underpinning molecular mechanism of those toxic effects have been yet unknown. CASQ2 possesses 
several indentations of different sizes and charges in its highly negative surface. The majority of these 
pockets were hydrophobic while a large central junctional pocket was surrounded by charged amino 
acids. Based on their sizes, it was difficult to focus on one particular pocket because many of them 
indeed were suitable for binding of pharmaceutical compounds. Based on the presence of plausible 
pockets, we hypothesized that there could be differential affinities between strong and moderate binder 
for a given site, and there could be possibilities for suitably sized compounds to occupy multiple sites 
simultaneously. We further tried to probe these hypotheses by solving the appropriate drug bound 
crystal structures without any success. However we were able to crystallize rabbit skeletal 
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calsequestrin (rCASQl, PDB ID: 3V1W) in complex with a MPD. The amino acid sequences of 
CASQl and CASQ2 are highly conserved among human, canine and rabbit. Therefore we were able to 
qualitatively exploit the MPD binding profile of rCASQl crystal structure for evaluating the molecular 
docking result of canine crystal structure. In addition we were able to correlate these findings with the 
competitive inhibition of pharmaceutical drug binding by CAP and EGC, which was studied via 
molecular docking and isothermal calorimetry. 

2.2.1. Prediction of Binding Pockets for Different Classes of Drugs — Molecular Docking 

The docking algorithm predicted that a large pocket located in the central junction, and a few other 
cavities located in the inter-domain junctions are the only sites accessible to drug molecules. We 
classified these sites as the SI, S2 and S3 binding pockets (Figure 3A-C). The SI pocket encloses 
junctional cavity formed between the three domains of CASQ, and is comprised of an inner narrow 
sub-cavity (Sl-IC; Figure 3F) and a larger outer cavity (Sl-OC; Figure 3A-C,F). These cavities are 
fonned by the amino acids Q^^^ l'^^ Y^^^ Q''', R''\ S''\ R''\ R''\ L''\ R''\ T''' (Sl-IC); 
and I^''^ E^^^ D^''^ F^™, R''\ R''\ E''\ D''' (Sl-OC). The S2 pockets (Figure 3A-C) are 
cavities present in the thioredoxin-like folds of each domain, and are comparatively smaller than the 
SI pocket. Although there are three S2 pockets, one located in the thioredoxin-like fold in each of the 
three domains, they displayed significant variations in size and surface charge due to variations in 
amino acid compositions within these pockets. The S2 site of domain III, located in the C-terminus, 
revealed the presence of a bigger hydrophobic pocket with long stretch of highly conserved amino 
acids, while that of domain I, located in the A^-terminus, was short and made up of only 8 amino acids. 
Of these, only the S2 sites of domain I and domain III were found to be preferred by the studied drugs. 
The S3 sites are the interfacial crevices formed between domains I, II and III (Figure 3A-C), but only 
one of them were predicted to host a few poor binding molecules; none were preferred by any strong 
or moderate binders. The S2 sites in domains I & III are formed by the amino acids: K^"^, S^"^, E^^^, F^^^, 
j^iOT^ Qi08^ pii8^ g243 ^^^^^ j (S2-D1); M^'^^ M''\ F''', V^", A^^^ F^^^ A^^^ 

g256^ j287^ j3289^ ^290^ p292^ p293^ ^294^ ^299^ g300^ -p302^ p303^ j305^ j3306^ ^^07^ p310^ q311^ j312^ q313^ y3U 

y3i5^ ^316^ y3i7^ j^3\9^ ^320^ j332i^ g322^ ^324 Q^domain III (S2-D3). The amino acid composition of the 
S3 site predicted to be accessed by poor binders is: D^^^ P^^^ E^^^ H^^^; S3-III/II: E^^^ H^^^ R^^', 
T''\ L'''; S3-IIW: E'\ F'"", D''\ G''\ E''', V"^ R^^^, D^^^ M^^^ E''' (S3-II/I). 

2.2.2. Strong Binders 

Molecular docking predicted that the strong binders, such as DAN and DOX, preferred to bind into 
the Sl-OC site in their top ranking poses, wherein their bulky aromatic rings were able to fit into the 
larger cavity and are stabilized by 7i-stacking. In addition, neighboring residues, H^^^, T^^^, and R^^^ 
established multiple hydrogen bonds, significantly strengthening their affinity for the SI site relative to 
the moderate and poor binders (Figure 3D,E). On the other hand, the strong binders, such as AMD and 
VPM, preferred the S2 site (Figure 3 A), while TDZ was predicted to bind into Sl-IC in its top-ranking 
conformation. The docking scores for those strong binders correlated well with (R = 0.80) the 
corresponding experimental ^d's determined by ITC studies [9]. The common binding mode observed 
in the top ranking conformations of DAN, DOX, and TFP (Figure 3F) suggested that they might have 
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similar electrostatic charge distributions and orientation of the steric bulk, which were in good 
agreement with the previously published quantum mechanics data calculated using minimal basis sets 
without any polarizable functions [8]. We also compared the charge distributions of the optimized 
structures of the predicted SI pocket binding drugs DAN, TDZ and TFP (strong binders) with DTZ 
(moderate binder) calculated using polarized basis sets, which is discussed later. A good overlap of top 
ranking conformations and charge distributions of DAN, TDZ, and TFP adds strength to the prediction 
that the Sl-OC is the most preferred pocket for known highly toxic pharmaceutical compounds. 

Figure 3. Binding-pockets predicted for the drug molecules. The different cavities and 
sub-cavities are highlighted by boxes with different colors and are individually labeled. 
(A) Strong binders; (B) Moderate binders; (C) Poor binders; (D) The predicted binding 
mode of doxorubicin (DOX) reveals its occupation at the larger Sl-OC, being strengthened 
by several H-bonds with the neighboring amino acids; (E) The predicted binding mode of 
diltiazem (DTZ) reveals a unique orientation in which it occupies both the Sl-IC (phenyl ring) 
and the Sl-OC; (F) Two boxes with dashed lines indicate different regions within the SI 
site. Sl-IC is indicated by the yellow line and Sl-OC is indicated the black line. These 
figures were generated using Open-Source PyMOL™ (vl.4). 




2.2.3. Moderate Binders 

The moderate binders AMT, DTZ, NTP, and PMT were predicted to preferentially bind into Sl-IC 
with a nice overlap of their aromatic rings in their top-ranking poses (Figure 3B). This is also 
supported by the fact that these molecules have very similar experimentally determined values [9]. 
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Visualization of docked results highlighted that the Sl-IC pocket was more preferred by these 
compounds than the Sl-OC due to their smaller aromatic ring sizes, compared to the bigger rings 
possessed by strong binders like DAN and DOX. Moreover, these molecules had fewer functional 
groups and were not suited to bind into polar Sl-OC cavity. DTZ was unique among the moderate 
binders in that its benzyl ring occupied the Sl-IC site and the thiazepine moiety occupied the Sl-OC 
(Figure 3E), indicating that when two loosely connected aromatic rings are present, one of them goes 
into the Sl-IC and the excess bulk is pushed into the bigger Sl-OC cavity. Among such two ring 
systems, the one with polar functional groups tend to go into the polar Sl-OC cavity, thus providing 
better opportunities for stabilization via H-bonded interactions. This hypothesis is supported by the 
fact that the thiazepine moiety of DTZ was predicted to form multiple hydrogen bonds with 1146, 
Q226 and T173, while the other moderate binders like AMT, NTP, and PMT showed no hydrogen 
bonding. The placement of the thiazepine moiety in Sl-OC resembled the common binding mode 
observed in strong binders, but, unlike DOX, DAN, and TFP, the DTZ it did not exclusively utilize the 
Sl-OC, explaining its ability to achieve affinity in the range of strong binders like DAN and DOX. 
Other moderate binders, IMP and NTP, were also predicted to occupy the Sl-IC, while AMT and PMT 
preferred either S2-DI or S2-DIII sites (Figure 3B). The reduced toxic side effects by these molecules 
suggest that stabilizing the S2 sites does not affect the biological function of CASQ as strongly as 
stabilizing the Sl-OC by DAN and DOX. 

2.2.4. Poor Binders 

Four of the five tested poor binders, EPH, IPT, MNZ, and MPL, preferred to bind either into one of 
the S3 pockets or S2-DII (Figure 3C), while HPD was predicted to bind to Sl-IC (Figure 3C). They 
did not reveal any preference towards the polar Sl-OC, which explains their poor binding nature due to 
lack of H-bonds. Nevertheless, the docking scores of those poor binders revealed higher scores in the 
range of moderate or strong binders and showed poor correlation with experimentally determined K^'s, 
which is possibly because of limitations in the docking method. It was previously shown at one 
instance that a basic amine group in some drug compounds did not formed any H-bonded interactions 
in the crystal structure, but molecular docking favorably top ranked the poses in which these groups 
formed H-bonded interactions with binding site residues [23]. Such circumstances lead to scoring bias 
and result in false positives with higher scores, as could have happened with our poor binding 
compound set. 

2.3. Quantum Mechanics Studies ofLigand Charge Distributions: Comparison between Strong 
Binders Predicted to Bind into the Sl-OC ofhCASQ2 

The low micro molar range [9,10] binding affinities of strong binding pharmaceutical drugs: DAN, 
DOX and TFP in CASQ2 have been reported and are well known for their cardiotoxic side effects. We 
wanted to know if the charge distribution pattern between the predicted Sl-OC binding drugs are 
similar and in accordance to the previously published data [8]. To accomplish this, we performed 
accurate quantum mechanical (QM) calculations for the strong binding drug molecules: DAN, TFP, 
TDZ; and a moderately binding molecule: DTZ in order to compare their optimized geometries and 
charge distribution properties (Figure 4). The QM calculations revealed that the strong binders TDZ 
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and TFP had planar phenothiazine moieties that were highly neutral in the core. The side chains of 
TDZ and TFP showed variations with some polar atomic substitutions. But, the surface of DAN 
revealed that the core anthracycline ring was neutral in accordance with the observations in TDZ and 
TFP, but highlighted plenty of terminal negatively charged functional groups. This infers that all these 
drug molecules involve bulky hydrophobic interactions for a tight binding, but the DAN could produce 
both hydrophobic and polar interactions. The structure of DTZ revealed difference in the out-of-plane 
aromatic rings, where the methoxy phenyl ring can rotate flexibly, as it happened in the optimized 
structure (Figure 4 bottom right). These results displayed the similarities in electrostatic charge 
distributions between them, which in turn explained the docking predictions that they all bind into an 
identical pocket. 

Figure 4. The iso-density surface maps for predicted hCASQ2Sl-0C binders using the 
density functional theory method employed in Gaussian 03 [24]. The red and blue areas 
correspond to electron-rich (negative) and electron-poor (positive) regions, respectively. 
The color bar located at the bottom represents the charge scale. 




Thioridazine Trifluoperazine 




Daunomycin Dlltlazem 



2.4. Amino Acid Conservation Analysis of Potential Small Molecule Binding Pockets in CASQ 

Amino acid sequence conservation profile of rCASQl was studied using the ConSurf [25-27] 
algorithm as employed in the web server. The target sequence (rCASQl) was used for multiple 
alignments with other homologous CASQs that were retrieved separately using the BLAST [24] 
algorithm. The conservation profiles of amino acids at each position were later quantified in the score 
range of 1 (no conservation) to 9 (fully conserved) by ConSurf This analysis revealed that at least 
60%-70% of the whole CASQ sequence is conserved among a wide range of species, while >90% of 
the pockets that were predicted to host the known sets of strong, moderate and poor binding drugs 
{i.e., SI and S2) are highly conserved (Figure 5C). Since these two sites represent the junctions 
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between the three domains of CASQ and can significantly affect the polymerization of CASQ. The 
overall conservation of amino acid sequences and both the SI-OC and SII pockets indicate that they are 
significant to the function of the protein and hence its normal function will be disturbed upon 
drug-association. 

Figure 5. Binding sites for Drug Molecules: Evidences from crystallography and 
predictions from molecular docking. (A) Surface representation of the rCASQl-MPD 
complex crystal structure. The three crystallographic 2-methyl-2,4-pentanediol (MPD) 
molecules are represented by ball and stick (green for carbon atoms and red for oxygen 
atoms). The corresponding electron density maps for MPD molecules were contoured at 
1.5 a level and are represented by sky blue mesh. The cavities accessible to MPDs as 
predicted by docking are represented by the blue balls; (B) Ribbon diagram representing 
the sites accessible for MPD molecules in the crystal structure of rCASQl (red and green) 
and for the predicted MDP molecules (blue and red). These figures were generated using 
Open-Source software PyMOL™ (vl.4); (C) Analysis of the sequence conservation of 
rCASQI. The evolution of each amino acid forming the potential binding sites are 
individually scored by the server ConSurf The average conservation scores are displayed 
at the far right hand side end, representing the average for each of the predicted binding site. 
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2.5. Structural Basis of 2-Methyl-2,4-pentanediol (MPD) Binding into rCASQl: Crystal Structure 
Reveals Simultaneous Multiple Binding into S2-DIII Pocket Predicted by Molecular Docking 

Several initial attempts to crystallize CASQl or CASQ2 in their drug-bound form did not succeed 
in the presence of drugs like DOX, DAN or TFP. However our attempt to co-crystallize rCASQl in 
complex with MPD, a molecule much smaller than the studied set of pharmaceutical drugs, generated a 
diffraction quality crystal (Figure 5A,B, PDB ID: 3V1W). The small and aliphatic nature of MPD 
presents it a high level of conformational flexibility and hence we expected it to simultaneously 
occupy a number of different surface pockets inrCASQl. The crystal structure clearly revealed the 
positions for three MPD molecules in agreement with our hypothesis of simultaneous multiple drug 
binding (Figure 5A,B). All three MPD molecules were found localized within the S2-DIII domain at 
the C-terminus region of CASQ that is highly hydrophobic in nature. This crystal structure provided an 
initial idea about possible small molecule binding pockets. 

2.6. Development of a Structure Model for Competitive Inhibition of Pharmaceutical Drugs by 
Caffeine (CAP) and Epigallocatechin (EGC) 

2.6.1. CAF/EGC Mediated Competitive Inhibition — Molecular Docking Studies 

We further developed a structure model to understand the CAF and/or EGC mediated antagonistic 
activities [8,22] against drug-induced muscular toxicities at the molecular level. To accomplish this, 
we docked the CAF, EGC and EGCG into CASQ by applying the same docking approach as the 
pharmaceutical drugs. The predicted binding sites for CAF and EGC from the top ranking binding 
poses were then compared with those of strong and moderate affinity drugs to check for any overlap. 
As expected, there were direct structural overlaps between CAF, EGC, and most of hCASQ2-affinity 
drugs in their binding sites (Figure 6A,C), suggesting a plausible mechanism for CAF/EGC-mediated 
competitive inhibition. The docking scores for CAF and EGC were — 3 and — 6 kcal-moF^ indicating 
a loose binding that correlated to poor and moderate binding range at S 1 and S2 cavities, respectively. 
The top scoring pose of CAF revealed binding into the SI pocket, while energy difference for it to 
bind into a S2 pocket was -1 kcal-mor\ However, EGC was predicted to bind into S2 cavity in its 
highest-ranking conformations and the energy difference was -0.1 kcal-moF^ for it to bind into the SI 
pocket. Such low penalties associated with occupying one or other pockets suggested possibilities for 
the simultaneous binding in more than one cavity for these natural compounds, in which case, all drug 
molecules that bind either to SI and S2 sites would be antagonized by them. The smaller number of 
aromatic rings and the hydrophobic nature seemed to allow CAF to fit snugly into the Sl-IC, but the 
size of its aromatic core extends out of the Sl-IC pocket. Thus, the docking results of CAF suggest a 
steric hindrance for strong binders like: AMD, DAN/DOX, and TFP to be accommodated in Sl-OC 
(red molecule in Figure 6A,C, Figure 7, and Supplemental file 2). The predicted CAF binding mode 
leads to: (i) the extension of bulky aromatic core into the Sl-OC cavity, and (ii) a high possibility of 
binding pocket rearrangement because of the observed tight fit of CAF core into the static crystal 
structure. Similarly, when the EGC bound into SI, it oriented in such a way that the trihydroxyphenyl 
ring occupied the Sl-IC and extended its fused ring into the Sl-OC in a manner highly similar to that 
of DTZ (Figure 6B), thus directly overlapping with the strong binding drug molecules. We further 
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addressed the possibilities for simultaneous multiple binding of CAF and EGC by checking their 
inhibitory activities that prevent the association of drugs that were predicted to bind inSl and S2 
cavities via ITC assays. 

Figure 6. The structural insights for competitive inhibition. (A) Overlap of the binding 
modes for caffeine (red)/gallocatechin (blue) and strong binders; daunorubicin (magenta), 
doxorubicin (brown), and trifluoperazine (yellow); (B) Overlap of the binding modes for 
diltiazem (green), epigallocatechin (blue), and epigallocatechin gallate; (C) Overlap of the 
binding modes for amlodipine (teal) and caffeine (red). These figures were generated using 
Open-Source PyMOL (vl.4). 




2.6.2. CAF/EGC Mediated Competitive Inhibition — Isothermal Calorimetry Assays 

We performed ITC assays to check the predicted competition between the natural compounds 
(CAF, EGC, and EGCG) and tested synthetic drugs. However, addition of EGC or EGCG to the buffer 
instantly precipitated the other components, limiting the inhibition studies to be carried for CAF only 
(Figure 7 and Supplemental files 2, 3). Among all of the studied drugs, the solubility of amlodipine 
was least affected by CAF, making it the most amenable drug for use in ITC competitive-inhibition 
assays. As shown in Figure 7, a series of titrations containing 50 |liM hCASQ2 solutions pre-incubated 
with 0, 0.5, 1.0, 2.5, and 5 mM of CAF showed a gradual reduction in the of amlodipine from 12 to 
33, 46, 66, and 144 |iM, supporting the computational prediction that CAF inhibits drug-binding to 
CASQ2. Although inhibitory effect of CAF to the strong and moderately binding drugs were 
established at a 5 mM CAF concentration, the concentration dependency could not be established for 
many other drug molecules due to unfavorable interactions with the solution that led to their 
precipitation. Data showing the inhibition by CAF for other CASQ2 affinity drugs are given in 
Supplemental file 2. Thus the ITC assays confirmed the docking prediction of competitive inhibition 
by showing significantly reduced affinity of all tested drug molecules to hCASQ2 in the presence of 
CAFITC is also suggestive of CAF binding into both S 1 and S2 cavities simultaneously, as a result of 
which, the drugs predicted to bind into SI and S2 cavities were successfully inhibited. 
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Figure 7. Comparison of the heat released upon amlodipine binding to hCASQ2 in 
the presence of 0 (pink), 0.5 (green), 1.0 (black), 2.5 (blue) and 5 mM (red) caffeine 
concentrations. The corresponding values for amlodipine are indicated. 
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2.7. Implications of Drug Binding on the CASQ Function 

We have previously shown that strong- and moderate-affinity CASQ-binding drugs significantly 
disrupt polymerization of CASQ (Figure 8) and consequently, its high capacity Ca -binding [9,18]. 
Our previous data with TFP indicate that Ca^^ is being released from SR microsomes upon its diffusion 
into SR [10,1 1] probably due to its association with CASQ and increased free Ca^^ concentration inside 
the SR above a threshold value of RyR channel opening, where similar sized hydrophobic control 
drugs showed no effect on Ca release [10,1 1]. Moreover, TFP was shown to significantly reduce the 
total SR microsomal Ca content in a concentration- and time-dependent manner [10,11]. As the 
predicted drug binding pockets, SI and S2, are located at the dimer interface of CASQ, we hypothesize 
that the association of drugs into either of these sites would significantly affect the interfacial 
geometry and consequently disrupt the polymerization of CASQ (Figure 8-B), which is a very crucial 
mechanism in the SR Ca^^ regulation. In addition, altered polymerization caused by drug binding 
could change CASQ's affinity to the Ca -release channel complex proteins such as RyR, triadin, 
and junctin, and, consequently, affect the opening probability of RyR2, although it is still highly 
controversial which oligomeric form of CASQ2 interacts with RyR2 [6,28,29]. Also substantial 
amount of polymorphisms have been known in CASQ2 sequence and its post-translational 
modifications [14,30]. This led to a general notion that certain individuals might be more sensitive to 
certain drugs than others in terms of muscle toxicity. It is tempting to speculate that some of those 
polymorphisms of CASQ could further increase the affinity of CASQ to certain drugs or their 
metabolites due to the altered characteristics of their binding pockets. 
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Figure 8. Schematic diagram of the polymerization of CASQ2. (A) The mechanism of 
Ca^^-sequestration via polymerization of the CASQ monomers; (B) The hypothetical 
inhibition mechanism of drug preventing polymerization of CASQ. The drugs are 
represented as big orange spheres; (C) In the presence of high concentrations of weak 
binders such as CAP and EGC (small yellow spheres), occupation by strong binders is 
being inhibited, while the weakly binding small molecules do not significantly affect 
CASQ2 polymerization. 
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3. Experimental Section 

3.1. Preparation of Native Rabbit Calsequestrin (rCASQl) 

Extraction of rCASQl from isolated SR microsomal was performed as previously described [10], 
but with an increase from 5 mM to 20 mM imidazole in buffer A and addition of 0.05% (w/v) NaNs. 
The SR microsomal extract was dialyzed into buffer B (20 mM MOPS, 0.5 M NaCl, 1 mM EGTA, 
0.05% NaNs, pH 7.2) and applied to phenyl-sepharose resin (GE Healthcare). After washing the resin 
extensively with buffer B, rCASQl was eluted using buffer C (20 mM MOPS, 0.5 M NaCl, 15 mM 
CaCl2, pH 7.2, 0.05% (w/v) NaNs), then buffer exchanged into buffer D (5 mM sodium phosphate, 
0.05% NaNs, pH 6.8) and applied to a CHT-IO-I hydroxyapatite column (BioRad, CA) equilibrated 
with the same buffer. rCASQl was then eluted from the column with a linear gradient from 60%-100%) 
buffer E (500 mM sodium phosphate, pH 6.8). Hydroxyapatite elution fractions containing rCASQl 
were pooled and exchanged into crystallization buffer (20 mM HEPES, 0.5 M NaCl, pH 7.0) and 
concentrated for crystallization. 

3.2. Crystallization of rCASQl and Structure Refinement 

Crystals of native rCASQl were grown at 4 "^C using the vapor diffusion method by mixing protein 
at a concentration of 10 mg-mL"^ in crystallization buffer with an equal volume of reservoir buffer 
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(0.1 M HEPES, 0.2 M NaCl, pH 7.0) containing 27.5% (v/v) 2-methyl-2,4-pentanediol (MPD) as the 
precipitating agent. Small crystals developed overnight and fully developed in one to two weeks. All 
crystals were of the orthorhombic space group C222i with one molecule in the asymmetric unit. 
Diffraction data were collected at the Advanced Light Source (Beamline 8.2.1) and reduced and scaled 
using HKL2000. Iterative adjustments and refinements were performed using COOT and PHENIX [31], 
respectively, using the deposited coordinates of native rCASQl (PDB ID: 3TRQ). The final 
coordinates have been deposited to the RCSB (PDB ID: 3V1 W), and were used for docking studies of 
MPD-binding sites. Complex crystallizations of CASQ from various species with several compounds 
such as daunomycin, doxorubicin, verapamil, diltiazem and trifluoperazine were attempted, but none 
diffracted well enough to be pursued. 

3.3. Iso-Density Surface Calculation 

Gas phase geometry optimizations using the B3LYP hybrid functional [32-34] and 6-31 + g(d) 
basis sets were performed in Gaussian 03 [35] in order to calculate the surface charge distribution for 
daunomycin (DAN), desipramine (DSP), diltiazem (DTZ), propafenone (PFN), thioridazine (TDZ), 
and trifluoperazine (TFP). The resultant electrostatic potentials were mapped on their van der Waals 
surfaces (Figure 1) using Gauss View 3 [36]. The surface charge distribution of drugs with 
similarly-charged surfaces were cross-compared to determine if they bound in a likewise manner in 
docking experiments. 

3.4. Isothermal Titration Calorimetry (ITC) 

ITC experiments were conducted using a previous protocol [9], with slight modifications. The 
protein and drug concentrations were kept at 50 |liM and 2 mM, respectively, while the CAF 
concentration was allowed to vary between 0 and 5 mM for different drugs. Compounds with poor 
aqueous solubility were dissolved first in DMSO and then diluted into ITC buffer with final 
concentrations of DMSO not exceeding 2%. Origin software (Origin Lab Corp, Northampton, MA, 
USA) was used to fit the data to a one-site binding model via non-linear least squares regression. 

3.5. Computational Molecular Docking 

All molecular docking simulations were performed using AutoDock 4.0 [37]. Ligand molecules 
were either retrieved from the NCBI PubChem web server [38] or generated with the Accelrys 
Discovery Studio Visualizer3.0 [39]. Any van der Waals clashes were cleared before final use. The 
ligand bond angles, bond distances, and dihedrals were treated as flexible in the docking process, while 
the protein was treated as a rigid receptor. Finally, these structures were converted into PDBQT format 
by AutoDock Tools [40], allowing only the polar hydrogen atoms to be present in docking. As the 
drug-binding site(s) in CASQ2 had not yet been established, blind docking for all drug molecules was 
performed by enclosing the whole protein inside a large grid box, to include all possible sites in the 
search. The atomic affinity grids were calculated using the AutoGrid module of AutoDock4.0. Default 
parameters were used, except those for torsional rotations, which were reduced from 50"^ to 25"^ in 
order to increase the exhaustiveness of the search. A total of 100 docking runs were performed for 
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each drug molecule using the Lamarckian Genetic Algorithm, resulting in 100 best conformations. 
Finally, the pool of predicted binding orientations were subjected to cluster analysis with a 2 A 
root-mean-square deviation cut-off, whereby the conformations within each cluster were ranked in 
order of increasing docking scores, which is a measure of the free energy of binding. 

4. Conclusions 

We report on the structural basis for drug-affmity of hCASQ2 by comparing molecular docking 
predictions with experimentally identified binding sites of multiple MPD ligands in the crystal 
structure of rCASQl. Our approach revealed possible binding pockets in hCASQ2 for known toxic 
pharmaceutical drugs. In addition, our molecular docking results indicated the possibility for 
simultaneous binding of multiple EGC and CAF molecules in more than one pocket, which was 
supported by the observation of multiple MPD-binding in the rCASQl crystal structure. The known 
drug molecules were sorted into three different groups based on their docking affinities to hCASQ2 as 
strong, moderate, and poor binders. Three different binding cavities: SI, S2, and S3 emerged as the 
binding sites for those compounds. Of these three sites, the polar outer cavity of SI (Sl-OC) was 
predicted to be the most preferred target for highly toxic drugs like DAN, DOX and TFP due to their 
bulky aromatic rings and opportunities for multiple hydrogen bond interactions with their polar 
terminal groups. However, these physico-chemical characteristics were relatively reduced in the 
moderately binding drugs. The known poor binding drugs had highly flexible side chains with several 
charged groups but were predicted to occupy the hydrophobic S2 or S3 sites due to smaller sizes of 
aromatic cores. Especially the most crucial difference seen in the poor binders were significant 
reduction in the size of their bulky aromatic rings. The CAF and EGC were predicted to bind loosely 
due to their smaller aromatic cores into the same binding pockets where the known highly toxic drugs 
were predicted to bind. Thus these natural compounds could compete with the drug molecules for the 
same binding pocket. These were consistent with our ITC data showing CAF mediated competitive 
inhibition of drugs binding to CASQ emphasizing the accuracy of our structure model. 
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